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Recent deep sequencing surveys of mammalian genomes have unexpectedly revealed pervasive and complex transcription and 
identified tens of thousands of RNA transcripts that do not code for proteins. These non-coding RNAs (ncRNAs) highlight the 
central role of RNA in gene regulation. ncRNAs are arbitrarily divided into two main groups: The first includes small RNAs, 
such as miRNAs, piRNAs, and endogenous siRNAs, that usually range from 20 to 30 nt, while the second group includes long 
non-coding RNAs (IncRNAs), which are typically more than 200 nt in length. These ncRNAs were initially thought to merely 
regulate gene expression at the post-transcriptional level, but recent studies have indicated that ncRNAs, especially IncRNAs, 
are extensively associated with diverse chromatin remodeling complexes and target them to specific genomic loci to alter DNA 
methylation or histone status. These findings suggest an emerging theme of ncRNAs in epigenetic regulation. In this review, 
we discuss the wide spectrum of ncRNAs in the regulation of DNA methylation and chromatin state, as well as the key ques- 
tions that needs to be investigated and acknowledging the elegant design of these intriguing macromolecules. 
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INTRODUCTION (Djebali et al., 2012; Lander et al., 2001; Venter et al., 


In 1956, Francis Crick described a “central dogma”, stating 
that genetic information flowed unidirectionally from DNA 
to RNA to protein (Crick, 1970). Around 20 years later, 
Paul and Duerksen surprisingly found that biochemically 
purified chromatin fractions contain more than twice as 
much RNA as DNA (Paul and Duerksen, 1975), raising the 
suspicion that RNA may be the primary organizer of chro- 
matin structure and epigenetic regulation. This underlying 
“dark matter” was extensively appreciated by transcriptome 
sequencing, which revealed that while ~1% of the mamma- 
lian genome carries protein-coding potential, more than 
75% of the genome is transcribed to produce a huge number 
of non-coding transcripts, including rRNAs, tRNAs, snR- 
NAs, miRNAs, siRNAs, piRNAs, and especially IncRNAs 
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2001). Although such pervasive and complex transcription 
is well accepted, the functionality of this vast number of 
transcripts needs to be further explored. Some think most 
ncRNAs are transcriptional noise or the by-products of 
transcription, whereas others claim they are functional 
(Hangauer et al., 2013; Jacquier, 2009; Wilusz et al., 2009). 
Moreover, RNA and chromatin are clearly intertwined to- 
gether, but mechanisms linking RNA to chromatin are still 
unclear. In this review, we first describe the transcriptional 
regulation of small ncRNAs and then focus on the epige- 
netic regulation of IncRNAs in diverse biological systems. 


ncRNA CLASSIFICATION: SIZE AND LOCATION 


There are three main types of classical RNA: messenger 
RNA (mRNA), transfer RNA (tRNA), and ribosomal RNA 
(rRNA). mRNA is able to be translated into protein, where- 
as rRNA and tRNA are housekeeping RNAs; they are high- 
ly abundant and play essential roles in translation process. 
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Non-coding RNAs are defined as functional molecules but 
without protein-coding ability. Based simply on their sizes, 
ncRNAs are arbitrarily divided into small and long RNAs 
(over 200 nt). Over the last two decades, several types of 
small non-coding RNAs, such as miRNAs, piRNAs, and 
endogenous siRNAs, have been discovered through genetic 
mapping. With the advent of deep-sequencing technology, 
an emerging small RNA world has been uncovered 
(Brosnan and Voinnet, 2009). These small ncRNAs play 
important roles in numerous cellular processes, such as cell 
proliferation, cell apoptosis, DNA damage, and induced 
pluripotency (He and Hannon, 2004). The abnormal expres- 
sion of small ncRNAs has been shown to cause tumors in 
many reports (Di Leva and Croce, 2010). 

The discovery of IncRNA occurred a bit earlier than that 
of miRNA, with H19 and Xist, involved in imprinting and 
dosage compensation, being the first two identified 
IncRNAs (Bartolomei et al., 1991; Brannan et al., 1990; 
Brown et al., 1991; Brown et al., 1992; Pachnis et al., 
1988). IncRNAs can originate from intronic, exonic, or in- 
tergenic regions (St Laurent et al., 2015). Based on their 
genomic location relative to the protein-coding genes, 
IncRNA has been further classified into four categories: (i) 
Antisense IncRNAs are transcribed in the opposite direction 
of protein-coding genes and have partial overlap with genes. 
(ii) Intronic IncRNAs are initiated in the introns of pro- 
tein-coding genes but without overlapping exons; (iii) bidi- 
rectional IncRNAs are defined by initiating in a divergent 
fashion from the promoter region of the nearest pro- 
tein-coding gene. (iv) Intergenic IncRNAs reside in inter- 
genic regions without overlapping any annotated pro- 
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tein-coding genes (Figure 1). Recent global run-on se- 
quencing in mammalian cells further identified a large 
number of IncRNAs derived from promoter and enhancer 
regions. Although their expression levels are very low and 
relatively unstable, active transcription of these RNAs is 
critical for optimal gene induction. Many of these enhancer 
IncRNAs belong to the intergenic IncRNA family, whereas 
promoter IncRNAs may belong to either the antisense 
IncRNA or bidirectional IncRNA classes (Rinn and Chang, 
2012). 


SMALL ncRNAs IN EPIGENETIC REGULATION 


Small regulatory RNAs, including miRNAs and siRNAs, 
are 20-23 nt small RNAs that were initially thought to reg- 
ulate gene expression by slicing mRNA or inhibiting trans- 
lation (He and Hannon, 2004). Currently, there are around 
1,800 miRNAs that have been identified in the human ge- 
nome. The number of miRNAs is still increasing exponen- 
tially along with advancements in high-throughput se- 
quencing technology (Griffiths-Jones et al., 2006; Kozo- 
mara and Griffiths-Jones, 2014; Luo et al., 2015). miRNAs 
were initially thought to only function in 3’UTR (untrans- 
lated region) regions, though recent studies show that 
miRNA functional sites are highly diverse in the cell (Chi 
et al., 2009; Helwak et al., 2013; Xue et al., 2013). Alt- 
hough thousands of miRNAs have been annotated in mam- 
malian genomes, defining their functional target genes re- 
mains a major challenge in the field. Thus, an efficient and 
unbiased high-throughput method would be valuable in re- 
visiting the principal of miRNA-mRNA interactions and 
re-evaluating the functional mechanism of miRNA. 
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Figure 1 The genomic landscape of RNA classification. IncRNAs are defined by their relative proximity to the nearest protein-coding gene. Coding exons 
are represented as cyan boxes, while non-coding exons are shown in other colors. Schematic representation of IncRNAs organized bidirectionally (purple), 
antisense to coding genes (red), in the introns of coding genes (rose), in intergenic regions (blue) or associated with a promoter (green). 
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Although miRNAs have not been reported to be directly 
involved in epigenetic regulation in mammalian cells, sev- 
eral groups have found that aberrant expression of miRNAs 
can alter the global DNA or chromatin state by restricting 
the activity of a related chromatin remodeling enzyme 
(Benetti et al., 2008; Denis et al., 2011; Garzon et al., 2009; 
Yuan et al., 2011). Considering the extensive base pairing 
of miRNAs with these pervasive transcribed nascent RNAs, 
it is not surprising that the miRNA pathway may directly 
participate in epigenetic control of gene expression. Indeed, 
two recent reports demonstrated the DNA methylation 
changes mediated by the miRNA pathway in plants 
(Khraiwesh et al., 2010; Wu et al., 2010). In line with this 
observation, genome-wide mapping of Argonaute 2 
(AGO2), the key component of the RNA interference si- 
lencing complex (RISC), revealed that over 60% of 
Ago2-miRNA sites are located in intronic, promoter, or 
intergenic regions (Xue et al., 2013), implying an undocu- 
mented role of miRNAs in transcription or epigenetic regu- 
lation. 

The silencing of repetitive regions, such as centromeres, 
telomeres, and the mat locus, by siRNA-mediated RNA 
interference has been well characterized in Schizosaccha- 
romyces pombe. Several groups have demonstrated that 
depletion of the key components of the RISC complex, such 
as Dicer (dcr1), Argonaute (Ago/l), and RNA-dependent 
RNA polymerase homolog Rdp1, causes aberrant accumu- 
lation of long non-coding RNAs from centromeric repeats 
and corresponding loss of histone H3 lysine 9 methylation 
(H3K9me), thus impairing the function of the centromeres 
(Hall et al., 2002; Volpe et al., 2002). Such long non-coding 
transcripts are not detected in wild-type cells, as Dicer is 
able to actively process them into small RNAs. These small 
RNAs in turn guide RNA-induced transcriptional silencing 
(RITS) complexes containing Agol, Chpl, and Tas3 to 
centromeric regions through base pairing with nascent tran- 
scripts or complementary DNA sequences (Noma et al., 
2004; Verdel et al., 2004). Meanwhile, the siRNA-mediated 
targeting tethers Clr4 methyltransferase complexes to repet- 
itive regions to methylate histone H3 at lysine 9. H3K9me 
in turn recruits chromodomain-containing proteins, such as 
Chp1, Chp2, and Swi6, to initiate the spreading and estab- 
lishment of heterochromatin domains. Interestingly, RITS 
complexes are proven to be able to interact with another 
complex called the RDRC (RNA-dependent RNA polymer- 
ase complex), which contains Rdp1 and Dicer. Further re- 
cruitment of RDRC is able to process nascent transcripts 
into siRNAs, which initiate another round of loading and 
tethering and thus form a self-enforcing loop to efficiently 
set up heterochromatin structure at the centromere and re- 
petitive regions in the genome (Motamedi et al., 2004). 

Such siRNA-mediated epigenetic regulation is also ob- 
served in Arabidopsis (Chen et al., 2014). In which a spe- 
cific argonaute family protein, AGO4, is responsible for 
locus-specific siRNA accumulation as well as DNA and 
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histone methylation (Zilberman et al., 2003). This RNA- 
directed transcriptional silencing and histone methylation is 
also conserved in mammalian cells. Several groups have 
demonstrated that transfected siRNAs targeting promoter 
regions are able to induce DNA methylation and histone H3 
lysine 9 dimethylation, although this is dependent on a dif- 
ferent argonaute protein family member, Agol (Han et al., 
2007; Kim et al., 2006; Morris et al., 2004). Moreover, re- 
cent massively parallel sequencing of small RNAs from 
mouse oocytes and embryonic stem cells identified a larger 
number of endogenous siRNAs and piRNAs, most of which 
are derived from retrotransposons, bidirectional transcrip- 
tion, and antisense transcripts; however, whether these 
siRNAs and piRNAs are involved in heterochromatin for- 
mation in mouse cells is still unclear (Babiarz et al., 2008; 
Watanabe et al., 2008). In line with these reports, an RNA 
component has been described to be essential for the 
maintenance of H3K9me-marked pericentric heterochroma- 
tin (Maison et al., 2002). Moreover, depletion of Dicer, a 
key component of the RNAi machinery, can cause signifi- 
cant defects in centromeric silencing in mouse embryonic 
stem cells and DT40 cell lines (Fukagawa et al., 2004; 
Kanellopoulou et al., 2005). Recent sequencing efforts re- 
vealed a strong correlation between H3K9 methylation and 
repetitive elements (Lippman et al., 2004), which account 
for two-thirds of human genome. These findings lead to an 
intriguing possibility that the RNAi machinery may have 
profound roles in epigenome maintenance and regulation. 
Further mechanistic studies of small ncRNAs involved in 
this process may help us to understand the inheritance of 
epigenetic regulation. 

Besides miRNAs and siRNAs, another small non-coding 
RNAs called piRNAs (piwi-interacting RNAs) have also 
been reported to play important roles in epigenetic regula- 
tion. In contrast to small RNA-mediated gene silencing, 
piRNAs were initially found to promote euchromatic his- 
tone modifications in Drosophila melanogaster (Yin and 
Lin, 2007). A recent genome-wide survey indicates that 
Piwi and piRNAs are extensively associated with chroma- 
tin, and that piRNAs are sufficient to recruit HPla and 
Su(var)3-9 to specific genomic loci to repress RNA poly- 
merase II transcription (Huang et al., 2013). This piR- 
NA-mediated repression has also been reported in other 
species, such as the Aplysia brain, in which a specific piR- 
NA is able to target Piwi protein complexes to the promoter 
region of the transcriptional repressor CREB2 and mediate 
its methylation at a conserved CpG island (Rajasethupathy 
et al., 2012). A similar mode of regulation may also exist in 
other somatic tissues, but this is still unclear. 


LncRNA IDENTIFICATION AND 
CLASSIFICATION 


The first identified IncRNA, H19, is one of the most abun- 
dant IncRNAs in developing mouse embryos; it is actively 
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transcribed from a maternal allele and reciprocally imprint- 
ed with insulin-like factor 2 (Igf2), which is paternally ex- 
pressed. Aberrant expression of H19 is closely linked to the 
genetic disorders such as Beckwith-Wiedemann syndrome 
and Silver-Russell syndrome (Bartolomei et al., 1991; 
Pachnis et al., 1988). Although it has been shown that H19 
plays important roles in cell proliferation and tumorigenesis, 
complete deletion of H19 has no observable phenotype in 
mice (Ripoche et al., 1997). Recently, two groups reported 
that H19 functions as the precursor of miR-675 or as a spe- 
cific sponge for let-7 family members, strongly arguing 
against its functionality as an IncRNA in imprinting balance 
(Cai and Cullen, 2007; Kallen et al., 2013). However, how 
H19 controls imprinting is still unclear. In 1991, another 
IncRNA, Xist was cloned from human X-inactivation re- 
gion. It originated from an inactivated X chromosome and 
was found to encode a 17-kb RNA that can coat the chro- 
mosome in cis. The essential role of Xist in X chromosome 
inactivation strongly suggests that a non-coding RNA itself 
could be a driving force in altering DNA and chromatin 
states (Lee, 2012). 

Around 15 years later, along with a consortium effort in 
epigenome mapping, ~1,600 IncRNAs were identified in the 
mouse genome by checking Pol II-transcribed genes with a 
signature of H3K4me3 and H3K36me3 chromatin features 
(K4-K36 domains) (Guttman et al., 2009). Now, the total 
membership of IncRNAs exceeds 200,000, although the 
number is still debated. This large number of IncRNAs im- 
mediately elicited questions about their functionality (Xie 
et al., 2014). However, IncRNAs show more evolutionarily 
conserved sequences than intronic and untranscribed inter- 
genic regions, suggesting that most of them are biologically 
functional in mammals (Bu et al., 2015). Moreover, Inc- 
RNA expression patterns are strongly associated with di- 
verse key biological processes, such as DNA damage, im- 
mune responses, and induced pluripotency, again indicating 
their potential functionality (Guttman et al., 2009; Song 
et al., 2014). In terms of the physical localization of 
IncRNAs relative to protein-coding genes, a subset of 
IncRNAs is predominantly localized in the nucleus, while 
the majority exclusively resides in the cytosol. Nuclear 
IncRNAs play important roles in epigenetic regulation, 
whereas cytoplasmic IncRNAs have been shown to regulate 
mRNA stability and translation. 


IncRNAs IN EPIGENETIC CONTROL 


Though a large number of IncRNAs have been identified, 
the mechanistic and functional roles for most IncRNAs are 
still unknown. Genetic studies have revealed that several 
IncRNAs, such as H19, Xist, Kcnqlotl, AIR, HOTAIR, and 
ANRIL, are associated with imprinting and heterochromatin 
formation. The common lessons from these IncRNAs are 
that they facilitate the recruitment of chromatin modifying 
enzymes to specific genomic loci to change the chromatin 
or DNA state (Rinn and Chang, 2012). Such actions are 
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employed by most IncRNAs. It has been reported that 
around 40% of IncRNAs are directly associated with diverse 
chromatin-modifying complexes, including PRC2, Co- 
REST, and SMCX (Khalil et al., 2009). The extensive 
number of associations between IncRNAs and chromatin 
modifying enzymes suggests in general that IncRNAs di- 
rectly guide chromatin-remodeling proteins to target loci 
and function as scaffolds for the recruitment of transcription 
factors to activate or repress gene expression. 

IncRNAs can function either in cis or in trans. For ex- 
ample, lincRNA-P21, a direct transcriptional target of P53, 
is activated by P53 upon DNA damage; it functions as a 
transcriptional repressor by recruiting the RNA-binding 
protein hnRNP K to the promoter region of P53 target 
genes. The detailed mechanisms of how lincRNA-P21 tar- 
gets hnRNP K to P53 target genes and how the whole com- 
plex inhibits transcription are still unclear (Huarte et al., 
2010). However, a recent study indicated that lincRNA-P21 
represses somatic reprogramming by further recruiting the 
H3K9 methyltransferase SETDB1 and the DNA methyl- 
transferase DNMT1 to the promoter region of pluripotent 
genes with hnRNP K, thus sustaining histone and DNA 
methylation levels and repressing pluripotent gene expres- 
sion (Bao et al., 2015). This hnRNP K-dependent recruiting 
mechanism may also be applicable to P53-dependent 
transcriptional repression (Figure 2). Like lincRNA-P21, 
HOTAIR is another well-characterized IncRNA involved in 
epigenetic regulation in a trans manner. It can simultane- 
ously bind two histone-modifying complexes, PRC2 and 
LSD1, via two highly structured RNA ends (Figure 2). This 
combination is coordinated with simultaneous H3K27 
methylation and H3K4 demethylation, thus ensuring the 
stable silencing of the HOX D locus (Tsai et al., 2010). As a 
strong indicator of breast cancer, elevated expression of 
HOTAIR is also coordinated with re-targeting of PRC2 
complexes to the promoter region of metastasis suppressor 
genes to silence their transcription by altering H3K27me 
levels (Gupta et al., 2010). Unlike HOTAIR, Xist and 
Kenqlotl mainly function in cis; they both interact with 
histone methyltransferase PRC2 to elicit locus-specific 
changes in histone H3K27 methylation (Rinn and Chang, 
2012). 

A recent genome-wide loss-of-function screening of 
IncRNAs in embryonic stem cells revealed that most 
IncRNAs function in trans, though this is still debated 
(Guttman et al., 2011). Several studies have indicated that 
cis or trans regulation is dependent on the relative expres- 
sion level of a specific IncRNA in a specific cellular con- 
text. For example, Xist is a well-characterized IncRNA that 
functions specifically in cis to regulate dosage compensa- 
tion; however, ectopic expression of Xist also induces de 
novo heterochromatin formation outside the X chromosome 
(Jeon and Lee, 2011). A similar phenomenon is also ob- 
served for HOTAIR (Gupta et al., 2010). Together, these 
observations suggest that how an IncRNA functions is also 
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Figure 2 Examples of functional mechanisms for several long non-coding RNAs. IncRNAs can function either in cis or in trans; for example, Kcnqlotl 
and Xist mainly function in cis to regulate chromatin remodeling. In contrast, HOTAIR and lincRNA-P21 primarily function in trans to repress gene expres- 
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strictly dependent on its relative concentration and stability. 
Once the concentration of an IncRNA reaches a sufficient 
level, it may acquire the ability to diffuse to ectopic sites to 
elicit global gene expression changes. 


THE EMERGING THEME OF AN 
RNA-CHROMATIN NETWORK 


As RNA is a critical component for chromatin organization 
and RNA is twice as prevalent as DNA in purified chroma- 
tin fractions, there is no doubt that RNA-chromatin net- 
works should the primary determinants of epigenetic control 
(Maison et al., 2002; Paul and Duerksen, 1975). In the fol- 
lowing sections, we will focus on the IncRNA Xist to illus- 
trate how these intriguing IncRNAs are designed to orches- 
trate crosstalk between the RNA and chromatin worlds. Xist 
is an X-inactivated specific transcript originating from the 
X-inactivation center (Xic). To balance X-chromosome 
gene expression between males and females, Xist is tran- 
scribed from the X chromosome to be silenced and coats the 
entire locus in cis (Figure 3). Once Xist is transcriptionally 
activated, it will recruit PRC2 complexes via the RepA el- 
ement. Meanwhile, the RepC element in Xist specifically 
binds the YY1 transcription factor, which in turn localizes 
Xist-PRC2 complexes to the nucleation center of the X 
chromosome and initiates its inactivation (Lee, 2012). 


Several models were proposed over the past two decades 
to describe how Xist-PRC2 complexes spread along the X 
chromosome for the de novo establishment of H3K27me3 
domains. To map Xist loading sites globally, two groups 
employed an antisense oligo-mediated affinity purification 
approach to identify Xist loading positions during X- 
chromosome inactivation. They revealed that Xist-PRC2 
complexes first target several gene-rich regions of the X 
chromosome to shut down transcription first and then 
spread to gene-poor regions locally (Simon et al., 2013). 
Importantly, jumping of the Xist-PRC2 complex from the 
nucleation center to distal loading sites seems to be de- 
pendent on the 3D genome architecture, as most of the 
loading sites are physically proximal to the nucleation cen- 
ter (Engreitz et al., 2013). Moreover, all the initial loading 
sites of the Xist-PRC2 complexes are located in transcrip- 
tionally active regions, indicating that transcribed nascent 
RNAs may be the initial signal for Xist-PRC2 targeting. 
The remaining question is how the Xist-PRC2 complexes 
are targeted to these nascent transcripts. 

Two recent reports provided initial clues for how 
Xist-PRC2 is loaded and how transcription is silenced at 
loading sites. Antisense affinity purification-mediated quan- 
titative mass spectrometry revealed several proteins that 
directly interact with Xist, many of which are RNA-binding 
proteins (Chu et al., 2015; McHugh et al., 2015). One of 
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Figure 3 The action model of IncRNA Xist in X-chromosome inactivation. Xist is transcribed from the X-inactivation center located on the X chromosome 
to be silenced. After activation, it will first recruit PRC2 complexes through the RepA sequence, and then transcription factor YY1 will bridge the complex 
onto the chromatin to methylate the future Xi. Later on, the Xist-PRC2/SHARP/HDAC3 complexes will jump from the X-inactivation center to initial load- 
ing sites, which were recently found to be gene-rich regions with high transcription activity. SHARP/SMRT/HDAC3 complexes will shut down transcription 
at initial loading sites and elicit histone H3 deacetylation. The Xist-PRC2 complex then will spread to gene-poor regions and methylate the entire chromo- 


some. 


these, called SHARP, is able to directly interact with the 
SMRT co-repressor and HDAC3, both of which are essen- 
tial for Pol II transcriptional silencing of the X chromosome 
(McHugh et al., 2015). Moreover, HDAC3 is able to initiate 
global histone H3 deacetylation in Xist territories, thus ex- 
plaining the long-standing observation of hypoacetylation 
on the entire silenced X chromosome, which appears very 
early upon initiation of X-chromosome inactivation. To- 
gether, these results raise the intriguing hypothesis that 
RNA-binding proteins in Xist complexes may provide an 
RNA code for the initial loading and targeting of Xist- 
PRC2. 


CONCLUSIONS AND PROSPECTS 


The ability of RNA-binding proteins (RBPs) to bind to non- 
coding regulatory RNAs and target them to specific nascent 
transcripts provides an ideal bridge between the RNA and 
chromatin worlds. A unique property of RBPs is that they 
can recognize specific RNA motifs, called RNA codes, in 
nascent RNAs. This may partially explain why a basal level 
of transcription is critical for IncRNA complex targeting. As 
IncRNAs tend to associate with several RNA-binding pro- 
teins, each of which has a preferred RNA sequence to bind, 


an array of different RBPs may acquire the ability to bind a 
long stretch of RNA sequences containing all of the corre- 
sponding RNA codes. This is likely the primary determinant 
that tethers IncRNA complexes to specific genomic loci 
(Figure 4). 


IncRNA 


Figure 4 Combinatorial RNA codes may determine the specificity of 
IncRNA/remodeling complex targeting. Each RNA-binding protein (RBP) 
has a specific RNA-binding code, which tends to be a 4-8-nt RNA se- 
quence (Ray et al., 2013). Several RNA codes in a specific arrangement 
may form an RNA code array. The corresponding RNA-binding proteins in 
IncRNA complexes may recognize the array containing RNA codes for all 
RNA-binding proteins; thus, the RBPs and nascent RNA sequence together 
act as the primary determinants for targeting of IncRNA complexes to 
specific chromatin loci. 
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IncRNAs, especially IncRNAs, are emerging as molecu- 
lar scaffolds to organize nuclear compartments and macro- 
molecular machinery in order to execute sophisticated and 
precise regulation of diverse biological processes. Under- 
standing how IncRNA complexes are organized and target- 
ed to specific loci, as well as where these loci are located, 
will be helpful in establishing new paradigms for epigenetic 
regulation. 
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